Facet Formation in the Negative Quenched Kardar-Parisi-Zhang Equation 
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The quenched Kardar-Parisi-Zhang (QKPZ) equation with negative non-linear term shows a first 
order pinning-depinning (PD) transition as the driving force F is varied. We study the substrate- 
tilt dependence of the dynamic transition properties in f +1 dimensions. At the PD transition, the 
pinned surfaces form a facet with a characteristic slope s c as long as the substrate-tilt m is less than 
s c . When m < s c , the transition is discontinuous and the critical value of the driving force F c (m) 
is independent of m, while the transition is continuous and F c (m) increases with m when m > s c . 
We explain these features from a pinning mechanism involving a localized pinning center and the 
self-organized facet formation. 
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The pinning-depinning (PD) transition by an external 
driving force has been of much interest recently. Typi- 
cal examples are interface growth in porous media under 
external pressure j), dynamics of a domain wall under 
random fields dynamics of a charge density wave 

under an external field [|J, and vortex motion in super- 
conductors under external current |^,^|. In the PD tran- 
sition, there exists a critical value F c of the driving force 
F, such that when F < F c , interface (or charge, or vor- 
tex) is pinned by disorder, while for F > F c , it moves 
with a constant velocity v. The velocity v plays the role 
of order parameter in the PD transition, which typically 
behaves as 



(F - F c f 



(1) 



Recently, several stochastic models for the PD transi- 
tion of interface growth in disordered media have been 
introduced @||. It is believed that the models in 1+1 
dimensions are described by the quenched Kardar-Parisi- 
Zhang (QKPZ) equation for the surface height h, 

d t h = vdlh+^{d x hf + F + r)(x,h), (2) 

where v and A are constants and the noise rj depends on 
position x and height h with the properties of (r)(x, h)) = 
and (r)(x,h)r)(x',h')) = 2DS(x - x')S(h - h'). The 
QKPZ equation for A > exhibits a PD transition with 
~ 0.64. The surface at F c can be described by the 
directed percolation (DP) cluster spanned perpendicu- 
larly to the surface growth direction in 1+1 dimensions. 
The roughness exponent a of the interface is given as 
the ratio of the two correlation length exponents of the 
DP cluster, that is a = vx/v\\ ~ 0.63. The coefficient 
A of the non-linear term renormalizes to A r which scales 
as X r ~ (F - Fa)-* (0 - 0.64). This is obtained by 
measuring the substrate-tilt dependence of the velocity 



where m is the tilt of the initial substrate || . 

Origin of the non-linear term in the QKPZ equation 
is different from that of the thermal KPZ equation with 
noise r/(x,t) [Q. For the quenched case, the non-linear 
term is induced by the anisotropic nature of disordered 
media, while for the thermal case, it is induced by the 
lateral growth effect and is proportional to the velocity of 
the interface which vanishes at the threshold of the PD 
transition [[fl]]. In particular, Tang et al. [[l0| showed in 
the context of vortex dynamics that the effective pinning 
strength takes the form (Ah + s 2 A x ) 2 / 3 /(l + s 2 ) 2 / 3 where 



s = d x h is the local slope and A, 1 / 2 and A x ' ' A are the 
amplitudes of random forces in the h and x directions, 
respectively. When the medium is anisotropic, Aft 7^ A x 
in general and the effective pinning strength depends 
on the local slope and generates the non-linear term in 
Eq.(||) with A cx (Ah - A x ). Therefore when A h < A x , 
i.e. when the surface is driven along the easy direction, 
A is negative. 

The negative QKPZ equation, i.e. the QKPZ equation 
with negative A is first studied in Jl^ ] . In marked contrast 
to the A > case, the negative QKPZ equation describes 
a first order transition in the sense that the velocity v 
shows a discontinuous jump at the PD transition. This 
is shown in Fig.[j]. At the PD transition, the pinned sur- 
face takes the shape of a mountain with flat inclination as 
shown in Fig.||. To elucidate the transition mechanism of 
the negative QKPZ equation, we study in this work the 
substrate-tilt dependence of the dynamic properties. We 
performed the direct numerical integration of Eq.(S 
one dimension with the discretized version used in 

h(x, t+At) = h(x, t) + At{h(x -l,t) + h(x + 1, t) 

-2h(x, t) + ^(h(x + 1, t) - h(x - 1, t)f + F} 
8 

+ (Ai) 2/3 £(z, [h(x, t)]), (4) 



1/2 



v(F,m) - v(F,0) + ym 2 



(3) 



where [• ■ •] denotes the integer part, and £ is uniformly 
distributed in [~, |]. The prefactor (At) 2 / 3 of the 
noise term arises from approximately coarse-graining the 
noise r](x, h) during a time interval At. This is different 
from that used in Ref. 113], but does not alter physical 
property of the PD transition. We choose At = 0.01 and 
use the initial condition h(x, 0) = mx and the helicoidal 
boundary condition h(L + i, t) = h(i, t) + Lm. Numerical 
results are discussed below in detail. 

First, we measure the critical driving force F and the 
facet slope at m=0 (no tilt) for several values of A. Fq is 
numerically defined as the value of F at which no surface 
in the sample of upto 200 configurations is pinned until 
a large fixed time to ~ 10 6 At. Just below F c , at least a 
finite fraction of the ensemble is pinned at to showing the 
facet morphology. The facet slope is determined by sam- 
ple averages of 100 configurations. Measured values of Fq 
and s c are shown in Table 1 for A = —0.5, — 1, —1.5 and 
—2. Also shown are vq, the velocity at Fo, i.e. the veloc- 
ity discontinuity at the first order PD transition. Next, 
we examine the tilt-dependence of the surface growth. 
For small tilts, the first order nature of the PD transition 
and the morphology of the pinned surfaces at the tran- 
sition do not change. Moreover, within our numerical 
accuracies, both the critical driving force and the critical 
facet slope are independent of the tilt m. We find this 
is so as long as m is less than a critical tilt m c which 
is essentially the same as the facet slope s c . When m 
becomes larger than s c , the surfaces at the pinned state 
cannot form a facet with slope s c by the helicoidal bound- 
ary condition and instead align along the substrate. In 
this case, the critical driving force F c (m) has to increase 
with m to compensate for the tilt-force. The schematic 
phase diagram in the F — m plane is shown in Figj^. 
The positively moving phase is bounded by a horizontal 
line at F = Fq for m < m c and a smooth curve F c (m) 
for m > m c . The transition to the v < phase (the 
lower curve in FigJ^) is the continuous PD transition of 
directed percolation universality class. This region is bet- 
ter understood by rewriting the QKPZ equation in term 
of h' = —h as 

dtti = vd 2 x b! - ^(d x h'f -F + V (x, ti). (5) 

The above equation implies that the negative QKPZ 
equation with negative F is equivalent to the positive 
QKPZ equation with positive F and vice versa. The 
negative velocity phase in F > region of Figj3| can be 
understood as a result of the positive QKPZ equation 
driven by a negative force which is not strong enough to 
overcome the effect of the positively driving non-linear 
term. We confine our discussion to the upper curve of 
Fig.|| below. 

A typical data set for the velocity as a function of F 
and m in the positively moving phase is shown in Fig^, 
while FigEI shows a three dimensional plot of v versus 



F and m. At F — Fq and for m < m c , v exhibits a 
discontinuous jump by Sv(m). The lowest curve in Fig.^ 
shows this. Since Sv(m) vanishes at m—m c , the point 
(F=Fo,m—m c ) may be regarded as a tricritical point in 
analogy with the equilibrium critical phenomena. Fol- 
lowing we measure A r by fitting the curves v(F,m) 
to Eq.(0) for small m. Since it is insensitive to F, we 
define A r as that obtained from F — F ; 

5v(m) = v(Fq, m) = vq — ^—^-m 2 + • • ■ . (6) 

The fifth column of Table 1 shows A r . For F larger than 
but close to F , v(F, 0) is linear in F. This together with 
Eq.(^J) allows one to approximate v(F, m) as 

v(F,m)=v +A(F-F )-^-m 2 + --- (7) 

where A is a constant. Setting v = in Eq.(|7]), one 
obtains an approximate form for F c (m) as 

F c (m)~F + ^(m 2 -m 2 c ) (8) 

with m c = (2t;o/|A, |) 1 / 2 . However, for larger m, v de- 
viates from Eq.(|7]) and higher order terms in m become 
significant. To find F c (m) and m c more accurately, we 
fit the velocity curve for large m and estimate F c (m) 
independently from the condition of v(F c (m),m) = 0. 
This is justified since 9 = 1 generically for m 7^ |1C|| . In 
Fig|| we show the values of F c (m) obtained in this way. 
The relation F c (m c ) = F Q then gives an independent 
estimate of m c which is also shown in Table 1. Note 
that m c obtained in this way agrees well with the facet 
slope s c as mentioned above. The values of m c itself 
strongly depends on A and could be fitted to a power 
law as m c ~ |A|~ a37 as shown in Figj7[ Thus when 
A = 0, m c becomes infinite, which recovers the previ- 
ous result that the critical force F c (m) is independent 
of the substrate-tilt in the quenched Edwards- Wilkinson 
universality class Jl(]||. On the other hand, |A r | scales 
as |A| - 8 in Fig|7| and hence m c - |A|-°- 37 - |A r |-°- 46 . 
This is to be compared with the zeroth order expression 
m c ~ |A r |-°- 5 in Eq.(|). 

In order to understand the origin of the discontinuous 
PD transition for A < 0, we examine the noise distribu- 
tion on the perimeter sites of critically pinned surfaces. 
It is found that the noise strengths around the site at the 
valley in Fig.|| are relatively small, corresponding to large 
pinning strengths. Meanwhile, when A < 0, surfaces tend 
to form a flat facet to make the term (X/2)(d x h) 2 nega- 
tively large, leading to a large pinning strength. However, 
the curvature force V 2 h is almost zero on the hillside, 
but is positive at the valley. Thus the sites around the 
valley are pinned by balancing the driving forces due to 
the curvature and the external force F with the pinning 
forces due to the noise and the non-linear term. Once 
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the sites around the valley is pinned, the stable slope s c 
of facet is selected in a self-organized way: Suppose that 
a surface is formed with a facet of slope s < s c . Then 
the velocity of the sites on the hillside would be positive 
according to Fig.[| The pinning strength due to (A/2)s 2 
is not strong enough to resist the external force, and the 
sites on the hillside move upward, and tend to make the 
slope of the facet larger. When the slope becomes larger 
than s c , however, the non-linear term is too strong, and 
the sites on the hillside move downward reducing the 
slope. Thus the slope s c becomes stable, and surfaces 
form a facet with the slope s c as shown in Fig.^. Since 
the surface pinning is initiated at the sites around the 
valley, once the pinning site is broken by an external 
force which is slightly bigger than Fq, the surface grows 
with a finite velocity, leading to the first order transition. 
On the other hand, when A > 0, the non-linear term is 
positive, and enhances the external driving force. Thus 
surfaces do not form a facet. The surface pinning is 
mainly due to local noise strengths as can be pictured in 
the DP cluster. Thus the pinned sites are not localized 
but scattered all over the system, so that the PD transi- 
tion is continuous. 

It would be interesting to observe the tilt-dependence 
of the negative quenched KPZ equation from a stochas- 
tic model. The Sneppen A model ]l2| , |l4j is known to 
belong to the negative QKPZ universality class. In the 
Sneppen A model, the restricted solid-on-solid (RSOS) 
condition is imposed on the height difference of neighbor- 
ing columns. Because of this restriction, the surface tilt 
cannot be larger than a maximum value. If the maximum 
tilt is less than the critical slope s c (A) of the stochastic 
model, the tricritical behavior would not be observed 
in simulations because the coefficient A is not control- 
lable. Recently a PD transition similar to the negative 
QKPZ equation is found in a seemingly different system, 
the driven Frenkel-Kontorova model |L5[ , where the PD 
transition is also discontinuous and the displacement at 
the pinned state looks like FigJ|. However, the tricritical 
behavior has not been studied yet in this system. 
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In summary, we have investigated the tilt-dependent 
behavior of the negative quenched QKPZ equation. We 
observed and explained that there exists a character- 
istic substrate-tilt m c such that the PD transition is 
discontinuous (continuous) when the substrate-tilt rn 
is less (greater) than m c . The characteristic tilt m c is 
found to be the same as the facet slope s c of the criti- 
cally pinned surfaces which are formed in self-organized 
way and depends on the non-linear term coefficient A 
as m c ~ |A|~ - 37 . Moreover, the threshold force F c for 
the PD transition is independent of m for m < m c and 
increases with increasing m for m > m c . The effective 
non-linear term coefficient A r remains finite as F ap- 
proaches F c . 

We would like to thank M. Kardar for suggesting this 
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FIG. 1. The velocity versus force to show the first order 
transition. The data are obtained from a flat substrate with 
A = —0.5. The solid line is a guide to the eye. 
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FIG. 2. Temporally evolved surface configurations of the 
QKPZ equation with A = —0.5 in 1+1 dimensions at the 
pinning-depinning transition point. Each curve is the surface 
profile at constant time intervals. The facet slope is s c m 2.1. 




FIG. 5. Three dimensional plot of the velocity versus sub- 
strate-tilt m and external force F for the case of A = —0.5. 
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FIG. 6. The depinning critical F c (m) versus the sub- 
FIG. 3. The phase diagram in the F-m plane. strate-tilt m for several values of A. 




FIG. 4. The velocity versus the substrate-tilt m for sev- 
eral external driving forces F for the case of A = —0.5. The 
curve at the bottom corresponds to the magnitude of the ve- 
locity jump at the critical force Fo, which becomes zero at a 
characteristic slope m c . 
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FIG. 7. Estimated numerical data for the critical tilt m c 
versus A. The broken line m c ~ |A| -0 ' 37 is a least square fit. 
Inset: The A-dependence of A r showing |A| ~ |A r | S . 



TABLE I. Measured values of the critical driving force Fo, 
the facet slope s c , the velocity discontinuity «o, the effective 
non-linear parameter A r , and the critical substrate-tilt m c . 



A 
-0.5 



Fo 
1.91 



2.08 



Vo 

1.20 



-0.62 



2.11 
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-1.0 1.99 
-1.5 2.06 
-2.0 2.09 



1.60 1.30 
1.40 1.35 
1.26 1.42 



-1.12 1.64 
-1.52 1.39 
-1.89 1.24 
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